Using unitary operations to preserve quantum states in the presence of relaxation 
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When a quantum system interacts with an external environment, it undergoes the loss of quantum 
correlation (decoherence) and the loss of energy (relaxation) and eventually all of the quantum 
information becomes classical. Here we show a general principle to use unitary operations to establish 
and preserve particular non-equilibrium states in arbitrary relaxing quantum systems. We elucidate 
these concepts with examples of state preservation in one-spin and two-spin entangled systems. 
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Quantum systems undergo damping and decoherence 
when they interact with thermodynamic systems. Pre- 
venting such non-unitary behavior is one of the biggest 
challenges facing the engineering of quantum information 
technology. 

Much work has been done to develop quantum error 
correction schemes to correct these non-unitary errors 
||. When the errors are below known thresholds 
these schemes can preserve quantum information for an 
arbitrarily long time. On the other hand, error correc- 
tion is likely to be difficult because the required fidelity 
is beyond current practice and a large supply of ancilla 
qubits is needed §§[§. 

When the dissipation modes are purely decoherence 
processes, Decoherence Free Subspaces that are unaf- 
fected by the decoherence can be relabeled and used for 
storing and processing quantum information (tJ Q . In 
this paper, we will look at the situation when the dissi- 
pation is damping or relaxation. Surprisingly in this case, 
unitary operations can be used to preserve a large sub- 
manifold of quantum states for an arbitrarily long time. 
We will outline the geometry of such situations and then 
detail two examples. The first example will illustrate the 
dynamics of one-spin systems. The second example will 
show how to preserve a pair of coupled spins in a partic- 
ular entangled state. 

First, consider a general open quantum system of di- 
mension TV with density matrix p. A general time evolu- 
tion of this system is given by the mapping 



P 1 



S t (p)=J2E(t)kpE(t)l 



(1) 



where the E k satisfy Yl k E k (t)^E k (t) = 1 for all t and 
are called Kraus operators p)). If our open system evo- 
lution is Markovian, and satisfies the algebraic property 
£t£ s = £t+s, then the evolution is called a quantum dy- 
namical semigroup. The dynamics are completely de- 
termined by the generator of this semigroup. The most 
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general differential equation for such a generator is given 
by the Lindblad Equation 



k=0 



L k ,pL{] + [L k p,L 



(2) 



where H is a Hermitian matrix and L k are a set of TV x N 
matrices If the L k are identically zero, then the 

Lindblad equation reduces to a Schrodinger equation and 
we call the matrix H the Hamiltonian part of the dy- 
namics. In turn, the set terms of the Lindblad equation 
involving the operators L k is called the dissipative part 
of the Lindblad equation. 

Denote the trace zero Hermitian matrices su(N) and 
choose a basis, F k , satisfying the orthogonality conditions 



T*(F 3 F k ) 
as a sum 



Every density matrix p can be written 



P : 



a + Efc=i r k F k 
N 



(3) 



If we change our representation and define the coherence 
vector, f — (r k ), then equation || can be written as a 
standard form ordinary differential equation 



Ar + Br 



(4) 



Here A corresponds to the Hamiltonian part of the Lind- 
blad equation while B and c correspond to the dissipative 
part 0. 

We will focus on a special type of quantum dynamical 
semigroups called relaxing semigroups. A semigroup is 
relaxing if for any initial state r(0), 



lim r(t) = Tf. 



(5) 



Such a situation occurs whenever the matrix A + B is 
invertible and the vector c is also nonzero. By setting 
dr/dt = 0, we see that ff = —{A + B) _1 c is the unique 
fixed point of the evolution. If the Lindblad equation 
has a unique fixed point, ff, then this fixed point is a 
global attractor. Indeed, in matrix form we find that the 
evolution is given in coherence vector form by 



f{t) = e ( A+B )*(fo-f/) + r> 



(6) 
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where ff is the global fixed point. 

Relaxing semigroups occur frequently in physical sys- 
tems where the equilibrium state of a system is known 
a priori. For example, in a liquid state nuclear mag- 
netic resonance experiment, the system will always re- 
turn to an equilibrium Boltzmann distribution which is 
solely a function of the applied magnetic field and the 
temperature. In what follows, we will show that when a 
quantum system evolves as a relaxing semigroup, unitary 
controllers can act to stabilize a variety of known states. 

Let us restrict attention to the following control sce- 
nario. Suppose a quantum system evolves as a relaxing 
semigroup, but that we can apply an arbitrary control- 
ling Hamiltonian, H c , to the system, but that we cannot 
adjust the dissipative terms. In the coherence vector rep- 
resentation we have 

r = A c f+Br + c (7) 

The Hamiltonian control cannot prevent relaxation, as 
the eigenvalues of the matrix A c + B will still have nega- 
tive real parts, but the following proposition shows that 
the controller shifts the fixed point of the relaxing semi- 
group. 

Proposition 1 Let B be the dissipative part of the Lind- 
blad equation for a relaxing semigroup. If B is diagonal- 
izable then for any Hamiltonian part A, the matrix A + B 
is invertible. 

Proof Let ( , ) be an inner product on . Then the 
skew symmetry of A implies that (r , Ar) = for all f. 
There exists a basis fj for B such that 

N 

(r,Bf)=J2(>>k+iu Jk )(rJ k ) 2 (8) 
fe=i 

where the < 0. The real part of the inner product is 
negative for any nonzero r. This in turn means that 

Re((r, (A + B)f)) < (9) 

for all r ^ which completes the proof. ■ 

If we apply a Hamiltonian A c then the state ff = 
[A c + B) -1 ? becomes the global attracting fixed point 
of our quantum system and A c is a stabilizing controller 
on our system. In particular, this means that when A c 
is applied, the system's steady state is ff, and and the 
system will flow to ff independent of the initial state. 
Hence the set 

C = {f = -(A + B) _1 c| A is a Hamiltonian} (10) 

can be made into fixed points of a relaxing semigroup 
using control Hamiltonians. Since these states will be 
stabilized by the dynamics, we will refer to C as the set 
of stabilizable states of our semigroup. 

We must note that the stabilizable states will in gen- 
eral be mixed states as the length of the vector ff will 



vary with the applied Hamiltonian. However, we will see 
that they can be useful for monitoring quantum systems 
and for preserving entanglement. The following theorem 
describes the geometry of the set C. In the proof, we will 
switch between the density matrix and coherence vector 
representations . 

Theorem 2 // the fixed point, p eq , of a relaxing semi- 
group has non- degenerate eigenvalues, then the set of sta- 
bilizable states is a simply connected N 2 — N manifold 
containing the fixed point of the process and having the 
maximally mixed state in its closure. 

Proof Let p eq be the fixed point of the quantum pro- 
cess with corresponding coherence vector f eq . Consider a 
small perturbation f = f eq + Sr. It is immediate to show 
that Af+b = ASf. 

For an infinitesimal time At, we have that p(At) = 
p + ASf At. If over this time, the eigenvalues of p(At) 
are the same as those of p(0), then there exists a unitary 
operator U with Up(t)W — p(0) and hence p is the fixed 
point of the process 

. . . exp(-iH c At)£ t exp(-iH c At)£ t . . . (11) 

which is generated by the Lindblad equation with H = 

Let \i/j n ) be an orthonormal eigenbasis for p eq with cor- 
responding eigenvalues p\ > ... > pn- We want to show 
that there is an TV 2 — N dimensional neighborhood of p eq 
where the eigenvalues are unchanged under such a small 
perturbation. Since p and £At(p) are perturbations of 
p eq , we can calculate the change in the eigenvalues 

Ap n = (lpn\Peg + Sp - £At(Peq + Sp)\lp n ) 



The set of matrices M G su(iV) such that (tj} n \M\tjj n ) = 
has dimension N 2 — N as it corresponds to those traceless 
Hermitian matrices with zeros on the diagonal. 

If p is in the set of stabilizable states and has a corre- 
sponding Hamiltonian H, then there is a corresponding 
Pn for pH. At the limit where p = oo, poc = 1/iV. 
Hence, t/N is a limit point of C. M 

The preceding argument sets an upper bound on the 
dimensionality of the space of stabilizable states. If the 
fixed point has degenerate eigenvalues or the set of con- 
troller Hamiltonians is restricted to a subspace of su(iV) 
then C will have smaller dimension. 

The utility of this formalism can be explored in a one- 
spin example. Consider the process of damping to the Z 
eigenket | f). 111 terms of the Bloch vector, the system 
will relax with time constant 71 = 1 jT\ along the Z-axis 
and decohere with time constant 72 = I/T2 > 71 /2 in 
the x-y plane. 

The Lindblad equation which generates such a semi- 
group is given by equation ^| with H = 0, K = 1, and 



L Q = ^I+ and L 1 = ^-^Z (13) 

where I + = (X + iY) / \[2 is the raising operator. In the 
form of equation || this amounts to the Bloch equations 
with ^4 = 0, 



B 




and C 




(14) 



Parametrizing the space of Hamiltonians in the Pauli 
basis gives the controller Hamiltonian 



A = 




(15) 



in terms of three parameter controls {u x ,u y ,u z } corre- 
sponding to rotations about the x,y, and z axis respec- 
tively. 

The fixed points of the Lindblad equation are given by 
the equation (A + B)r + c = which can can be solved 
to find the manifold C 



!=(*-£) a if) 



(16) 



C is an ellipsoid containing both the fixed point and the 
maximally mixed point as we proved earlier. Its minor 
axis is governed only by the ratio of 71 to 72 . 

We also find the appropriate open loop controllers to 
reach the state (x, y, z) on C are 



u x — i~2 — and u v = t\ — 

z z 



(17) 



To ground this example in practice, let us describe 
how it is readily applied applied to pulsed NMR. We 
can asymptotically reach the desired steady state by ap- 
plying the unitary operations exp(— i(u x X + u y Y)At) at 
a repetition rate of At. When T2 (the transverse relax- 
ation time I/72) is comparable with T\ (the longitudi- 
nal relaxation time I/71), the steady state component 
of the Bloch vector in the x-y plane can be asymptot- 
ically close to |. In the language of NMR, the steady 
state magnetization is equal to half of the peak magne- 
tization from a ir/2 pulse when T\ = T2. On the other 
hand, when T 2 ^ ?i, the steady-state magnetization ap- 
proaches zero. Ernst and Anderson jL2] and Freeman 
|fl3fl derived these steady state from the Bloch equations, 
and our current formalism includes their results as a spe- 
cial case. Furthermore, by varying the pulse width, and 
in turn the steady state, we have experimentally demon- 
strated control over the NMR magnetization vector over 
times much larger than Tl as shown in figure |l|. 

Investigations into steady-state NMR on multiple spin 
systems has been less broadly investigated. It has proved 
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FIG. 1: Experimental data tracing out the letters "ML" in 
the transverse magnetization of an NMR spin system. We 
prepared a Copper Sulfate sample in water following the pre- 
scription in [fL4| to create time constants Ti = T2 = 39 ± 5ms. 
The signal was measured using a Varian 500 MHz NMR spec- 
trometer. 



successful for the specific case of studying spin-lattice re- 



laxation in dipolar solids 
for multiple spins has not 



|15||lqj, but a general theory 
seen established. Our results 
readily extend to higher dimensional quantum systems, 
but parametrizing the set of stabilizable states and their 
corresponding controllers becomes much more difficult 
as the number of variables in the coherence vector scales 
quadratically with the number of levels. Nonetheless, we 
will demonstrate techniques for dealing with such larger 
systems and describe a particular example of using lo- 
cal controllers and an entangling operation to preserve a 
highly entangled state. 

The entanglement of a pure state, p = \tp){4>\ of two 
spin half particles is defined as 



E(p) = Tr(pi logpi) = Tr(p 2 logp 2 ) 



(18) 



where pj denotes the partial trace over the Hilbert space 
of the jth spin. The quantity E takes values between 
and 1 and provides an information-like measure of the 
entanglement between a pair of spins. Correspondingly, 
an "ebit" is a unit of entanglement. One EPR pair has 
one ebit of entanglement. For a mixed state, we can 
define the entanglement of formation to be the minimum 
amount of entanglement required to create this mixed 
states from pure states H . Precisely, 



E(p) = mmJ2Pi E (^) 



(19) 



where the minimum is taken over all combinations of pj 
and \ipj) which yield p — J2jPj\' l l ; j)(' t f } j\- Wootters found 
a functional form for this quantity which involves extract- 
ing the eigenvalues of an algebraic function of the density 
matrix [17j. Using this metric, we show how to construct 
a stabilizable state with an entanglement 0.355. 

To simplify the equations used to solve for the fixed 
point, we will restrict our attention to a simple model. 
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Consider a two-spin system where both spins undergo 
damping to the spin-up state identically and indepen- 
dently with constants 71 = 72/2. This corresponds to 
dissipative operators 



L o = \fl I+ ® 1 and Li = ^yl ® 7 H 



(20) 



Assume the spins are coupled via the Hamiltonian H = 
JZ\Z 2 . The fixed point of this evolution is the state 

l^o) = I TT). 

Allow for only local Hamiltonians to be applied. Then 
the admissible Hamiltonians can parametrized as 



u Xl X 1 +Uy 1 Y 1 +u Zl Z 1 +u X2 X 2 +Uy 2 Y+u Z2 Z 2 (21) 



The coherence vector is given by rjk = Tr((ij ® Cfcp) 
where a 3 = {1, X, Y, Z} for j = 0, 1, 2, 3. Putting this all 
together, we get the fifteen equations 







(22) 



with the coefficients of A and B readily solved for by 
algebra. 

Consider the Hamiltonian 

H c = ^-X x - JZ X + ^-X 2 - JZ 2 (23) 

and let |^) = | TT) and \fo) = (| U) + I it»//2. 
Inverting the system in equation ^2] and taking the limit 
as J approaches infinity, yields the fixed state 



(24) 



which indeed has the entanglement of formation of 0.355. 

The rate at which the fixed point approaches infinity is 
plotted in figure Even for relatively small ratios, J/7, 
this state is close to p e . In the context of quantum com- 
putation, this procedure could be used to make a "well" 
of entanglement. Spins that are coupled locally can be 



used to store a known entangled state and then their state 
can be swapped into another system which can process 
the entanglement for communication or computation. 

We have shown a method for analyzing relaxing semi- 
groups and have also shown that by applying control 
Hamiltonians the fixed points of these systems can be 
shifted. We have further demonstrated how to apply 
these techniques to preserve known quantum states for 
arbitrarily long times without the requirements of redun- 
dancy or error thresholds from quantum error correction. 

A framework for labeling and exploring the space of 
stabilizable states in higher dimensional systems remains 
to be determined. Already for two-spins there is no intu- 
itive description of the manifold of stabilizable states and 
we have only demonstrated one example of a state which 
can be stabilized. Combining this higher dimensional la- 
beling with a prescription for using the stabilizable states 
in a coherent fashion for quantum information process- 
ing could provide a new method for protecting quantum 
computers from thermodynamic errors. 
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FIG. 2: The entanglement of formation of the fixed point 
under the Hamiltonian of equation 03. 
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